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Abstract. The basic idea of an extremely fast convergent iterative method, 
the Forth-and-Back Imphcit Lambda Iteration (FBILI), is briefly described and 
the applications of the method to various RT problems are listed and discussed. 



1. Introduction 

Radiative transfer (RT) is the underlying physical phenomenon in many astro- 
physical problems and among the most difficult to deal with. The main difficulty 
arises from the non-local and in general non-linear coupling of the radiation field 
and the state of the gas. 

In the problem of NLTE line formation in a given atmospheric model, the 
internal state of the gas depends, via radiative transitions, on the radiation field 
intensity, which in turn depends, via the RT process, on the state of the gas 
over a wide range of distant points. Mathematically, the non-local coupling 
is performed by the simultaneous solution of the corresponding RT equation 
describing the dependence of the mean intensity on the source function J = 
A [5] by means of the so-called A operator and the statistical equilibrium (SE) 
equations defining the source function in terms of the mean intensity of the 
radiation field S = S{J). 

In the well-known two-level-atom line formation problem the non-local cou- 
pling is linear and the problem can be easily solved by using either direct or 
iterative methods. In a more general multilevel case, RT problem is non-linear 
and, therefore, an iterative method is required. 

The most straightforward iterative procedure, the so-called A iteration, 
solves radiative transfer and statistical equilibrium equations in turn. However, 
in most cases of interest (in scattering dominated media of large optical thick- 
ness) the rate of convergence of this simple procedure is infinitely slow. 

A broad class of ALI (Approximate (or Accelerated) Lambda Iteration) 
methods, currently in use, is based on the use of certain physical or computa- 
tional approximations of the A operator within an iterative procedure. These 
methods usually employ some free parameter controlling the convergence and 
almost always need additional acceleration by some mathematical techniques 
(Ng accelerati on, s uccess ive over-r elaxation method, etc.) to achieve high con- 
verg ence rate ([Hubenv 20031: [Atan ackovic-Vukmanovic 2 004) . 

lAtanackovic-Vukmanovic. Crivellari. &: Simonneaui (jl997l ) developed the Forth- 
and-Back Implicit Lambda Iteration (FBILI) - a simple, accurate and extremely 
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fast convergent method to solve NLTE RT problems. FBILI dramatically accel- 
erates the convergence of the classical A iteration while retaining its straight- 
forwardness. In this paper the basic idea of the method is briefly explained and 
the applications to various RT problems are shown and discussed. 

2. The Idea of FBILI Method 

In order to demonstrate the basis of the FBILI method wc shall consider the well- 
known case of the two-lcvcl atom line formation (with complete redistribution 
and no overlapping continuum) in a plane-parallel and static atmosphere. Under 
these assumptions, the specific intensity of the radiation field Ixij,{t) is described 
by the RT equation of the form: 

fi^ = Mlcc,ir)-S{T)] , (1) 

where r is the mean optical depth, x is the frequency displacement from the 
line center in Doppler width units, /x is the cosine of the angle between the 
photon's direction and the outward normal and ipx is the absorption-line profile, 
normalized to unity. The line source function (SE equation for a two-level atom) 
has the following form: 

S{T) = eB{T) + il-e)J^{T) , (2) 

where e is the standard NLTE parameter representing the branching ratio be- 
tween the thermal (LTE) contribution B{t) and the scattering term 

J<p{r) = J Jx'Pxdx = \J^diJ, J dxipxlxi^ir) , (3) 

which accounts for the angle and frequency coupling of the specific intensities 
at the given depth point r. 

The basic idea of the forth-and-back implicit A iteration in the solution of 
the problem is as follows. First, as suggested by the existence of two separate 
boundary conditions, the FBILI uses a separate description of the propagation 
of the in-going intensities of the radiation field /^^^(t) with initial conditions at 
the surface (r = 0) and of the out-going intensities Ix^{t) with initial conditions 
at the bottom of the atmosphere (r = tn)- Second, although the values of the 
radiation field are unknown, its propagation can be easily represented by using 
the integral form of the RT equation and assuming polynomial (e.g. piecewise 
quadratic) representation of the source function between two successive depth 
points. Thus, for each depth point one can write linear relations for the specific 
intensities as functions of the unknown values of the source function and of its 
derivative. Following the idea of iteration factors, it is the iterative computation 
of the coefficients of these implicit relations (implicit, as the source function is a 
priori unknown), rather than that of the unknown functions themselves, which 
greatly accelerates the convergence of the direct iterative scheme. 

In the first part of each iteration (forward process), proceeding from the 
upper boundary condition, using the integral form of the RT equation for the 
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in-going intensities for each layer {ti_i,ti): 

where At; = ti — t;_i, and assuming paraboUc behavior for the source function 
S'{ti-i) = 2{S{ti) — S{ti_i))/Ati — S"(r/), one can write the hnear local implicit 
relation 

I^^iri) = + b-^]S{Ti) + c-^S'{Ti) , (5) 

representing the values of the in-going intensities Ix^iTi) ^ given optical depth 
point Ti in terms of yet unknown values of the source function S{ti) and of its 
derivative 5"(r/). Here, = I~^(T;_i)e~^'^''^^/'^ + qS"{Ti-i) is computed with 
the old (known from the previous iteration) source function 5"°, whereas and 
c~^ depend only on the known optical distance Ar^. By integrating (5) over all 
frequencies and directions we obtain the linear relation 

J-{Ti)=b-S{Tl)+C~S\Ti) (6) 

representing implicitly the value of the in-going mean intensity. Thus, in the 
forward process, we differ from the classical A iteration that re-calculates J~ 
from the old (known) source function S°{t), in using the old source function 
to compute, at each optical depth point ti{1 = l,N), the coefficients bj^ and 
q~ of the linear relation (6). The coefficients are stored for further use in the 
backward process of computation of the new values of 5(r). Let us note here 
that the ratio in Eq. (5) of the non-local part of the in-going intensity to 
the current source function S°{ti) is actually the iteration factor. Since it is the 
only information that is carried from the previous iteration step, an extremely 
fast convergence is to be expected. 

The final aim is to derive, at each optical depth point ti, an implicit linear 
relation between the 'full' mean intensity and the source function 

MTi) = ai + biS{Ti) (7) 

that, together with equation (2), leads to the new source function. To obtain 
this, we need the coefficients of the corresponding relation for the out-going 
mean intensity. 

In the backward process we proceed from the bottom where I^^{tn), i.e. 
J^{tn) is known or more precisely, the coefficients of the implicit relation for 
the out-going specific intensities and, therefore, for the out-going mean intensity 

J^{Ti) = at + btS{n) , (8) 



are k nown at r = r^y. Eliminating S'{tn) from Eq. (6) (see lAtanackovic-Vukmanovic et al 



19971 ) we obtain the linear implicit relation (7), which together with SE equa- 



tion (2) leads to the new value of S{Tj\f). The new values of S{Tj\f), S'{Tj\f) and, 
hence, Ix^{tn) are then used to compute the coefficients of the linear relation 
(8) in the next upper layer. Together with the coefficients of Eq. (6) (stored in 
the forward process), we obtain the coefficients a/ and 5/ of Eq. (7). The com- 
putation of ai and bi and the solution of Eq. (7) together with SE equation (2) 
to get a new source function S{ti) are performed during the backward process 
layer by layer to the surface. The process is iterated to the convergence. 
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3. Applications of FBILI 

The accuracy and efficiency of the FBILI method have been checked in several 
RT problems. Here, we list some of the applications. 

3.1. Two-Level- Atom Line Formation 

The method was first developed for the two-level-atom line formation (with com- 
plete redistribution and no over lapping continuum) in a plane-parallel constant 

prop erty medium (for details see I Atanackovic- Vukmanovic 199ll . [Atanackovic- Vukmanovic et al 



119971 1. This case represents an ideal test for checking the numerical accuracy 
and stability of a new method. Namely, under such conditions the features of 
the solution depend only on the NLTE parameter e that is usually very small, 
so that the numerical errors can easily blur the solution. 

Thus, in order to test the stability of the method we solved the two-level 
atom problem with e = 10~^^. The results are gi ven in Fig. 1. They ar e 
compared with the exact discrete-ordinate solution ( Avrett fc Hummer] 1 19651 ) 



obtained by using the same discretization in optical depth (10 points per decade). 
The asymptotic value of the maximum relative error of the order of 0.3% is 
reached already within 9-14 iterations. Namely, only nine iterations are sufficient 
for the maximum relative correction between two successive iterations (for all 
depth points) to be less than 6 = 10~^ and 14 iterations for 6 = 10~^. Hence, 
we see that a negligible additional effort (of the iterative computation of the 
coefficients of the implicit relations, instead of the mean intensities themselves) 
with respect to the classical A iteration results in an extremely fast convergence 
(about 10 FBILI iterations compared to about 1/e classical A iterations, while 
one FBILI iteration takes only about 10% more CPU time than a classical A 
iteration) . 

In the literature the performances of various methods are usually given for 
e = 10~^. Thus, for this case, in Fig. 2 the convergence properties of the FBILI 
method are compared with those of t he ALI method s that use diagonal and 
3-diagonal approximate A* operators ( Hubenv 20031 ) . Excellent convergence 



properties of the FBILI are evident. Using FBILI yields convergence that is 
comparable to or even faster (for 6 > 10~^) than the 3-diagonal operator with 
Ng acceleration! 

3.2. Two-Level- Atom Line Formation: Partial Redistribution 

The FBILI method was applied to the case of the two-level atom line formation 
problem in which partial redistributi on is taken into account in the paper by 
Atanackov ic- Vukmanovic et al.l ( 19971 ). The results reproduced the well-known 



ones by iHummer (1969 ) . For the cases e = 10-"^ and e = lO"*, 13 and 15 
iterations, respectively, are enough to fulfill the criterion 5 = 10~^ for all the 
frequencies and all optical depths. 

3.3. Multilevel- Atom Line Formation 

The procedure for multilevel problem is the same as in the two-level-atom case. 
Starting with the known set of level populations, we repeat the entire forward 
process for each radiative transition z — ?> j. In the backward process, layer by 
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Figure 1. Evolution with iterations of the source function for a two- level 
atom with e — 10~^^ and B = 1. 




iteration 

Figure 2. Maximum relative change 6 of the source function as a function 
of the iteration number, for different variants of the ALI scheme (reproduced 
from Fig. 1 in Hubeny 2003) and for the FBILI method. In all cases, two- 
level atom without continuum is considered, with e = 10^'*, B ~ 1 and 4 
depth-points per decade of optical depth. 



layer, we compute the coefficients of the linear relation (7) for all the transi- 
tions, and replacing them in the SE equations we solve the latter for the new 
set of level populations. The test is performed by solving the same problem 
(three-level hydrogen a tom line formation in an isothermal atmosphere) as in 
Avrett &i Loeserl ( 19871 ). The solution with a maximum relative error below 3% 



is obtained in only nine iterations wi th the convergence criterion 5 = 10 (see 
Atanackovic-Vukmanovic et al.l 19971 ). 



3.4. Spherical Radiative Transfer: Monochromatic Scattering 

Th e generalization of the F B ILI m ethod to spherical geometry is performed 
by lAtanackovic-Vukmanovid ( 20031 ). Monochromatic scattering problem in a 



sp herical atmosphere is solved and the result s are compared with t hose given 
bv lGros. Crivellari. fc Simonneaul (|l997l ) and I Avrett k Loeseil (|l984l ). The rel- 
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ative difference of the solutions is about 1%. The solution is obtained already 
in the second iteration, whereas three iterations are required for the maximum 
relative correction to be less than 1%. 

3.5. Spherical Radiative Transfer: Line Formation 

In order to test the feasibility of the method when applied to the line forma- 
tion in spherically symmetric media, the test problem of the line transfer with 
background absorption , proposed by iMihalas. Kunasz. &: Hummerl ( 19751 ) and 



Avrett Sz Loeseil ( 19871 ). is solved. The solution is obtained in 15 iterations with 



an error less than 2%. 



4. Conclusions 

Forth-and-Back Implicit Lambda Iteration method is a simple, stable and ex- 
tremely fast convergent iterative method developed for the exact solution of 
NLTE RT problems. A negligible additional computational effort with respect 
to the classical A iteration results in an extremely fast convergence. No addi- 
tional acceleration is needed. The method is easy to apply. No matrix formalism 
is required so that the memory storage grows linearly with dimension of the prob- 
lem. Due to its great simplicity and considerable savings in computational time 
and memory storage FBILI seems to be a far-reaching tool to deal with more 
complex problems (multidimensional RT, RT in moving media, or when RT has 
to be coupled with other physical phenomena). 
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